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ABSTRACT 

In contrast with numerous studies on the physical and chemical structure of low- and high-mass protostars, much less is known about their 
intermediate-mass counterparts, a class of objects that could help to elucidate the mechanisms of star formation on both ends of the mass range. 
We present the first results from a rich HIFI spectral dataset on an intermediate-mass protostar, OMC2-FIR4, obtained in the CHESS (Chemical 
HErschel SurveyS of star forming regions) key programme. The more than 100 methanol lines detected between 554 and 961 GHz cover a range 
in upper level energy of 40 to 540 K. Our physical interpretation focusses on the hot core, but likely the cold envelope and shocked regions also 
play a role in reality, because an analysis of the line profiles suggests the presence of multiple emission components. An upper limit of 10 -6 is 
placed on the methanol abundance in the hot core, using a population diagram, large-scale source model and other considerations. This value 
is consistent with abundances previously seen in low-mass hot cores. Furthermore, the highest energy lines at the highest frequencies display 
asymmetric profiles, which may arise from infall around the hot core. 

Key words, stars: formation - ISM: abundances, kinematics and dynamics, molecules 
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1. Introduction 

Intermediate-mass, and therefore intermediate-luminosity proto- 
stars offer insights into the physical and chemical differences be- 
tween the formation of low- and high-mass stars, but questions 
about their chemistry and dominant gas-heating mechanisms re- 
main only partly answered. 

Deep in the interior of a protostellar core, energy is re- 
leased by a forming protostar. This energy heats the surround- 
ing gas by dust-mediated and UV photon heating, as well as 
through shocks caused by protostellar outflows. The respective 
roles of these mechanisms as a function of protostar luminos- 
ity and the effects of heating on protostar evolution, are un- 
der intense study (e .g. |Spaans et a"T1|1995[ |Doty et al.||2006| 
Brud erer et al.|[2009 1, with the role of shocks and UV photons 
in low- and intermediate-luminosity sources emphasized by re- 
cent Herschel Space Observatory results ( |Fich et al7||2010[ |van| 



The two main paths to this complexity are gas-phase and grain- 
surface reactions (e.g. Garrod et al. 2008). During the gradual 
warm-up of grains in a hot core, species formed on them in ear- 
lier evolutionary phases react and the products are later released 
into the gas phase. Observable chemical differences include a 
methanol-formaldehyde abundance ratio, which increases with 
decreasing protostar luminosity ( Cazaux et al.||200"3 Bottinelli 
|et al |2004||2"007l l. 

We present Herschel-HIFl sub-millimeter observations of 
methanol line emisson toward the intermediate-mass protostar 
OMC2-FIR^] attributing part of the emission to a hot core. 



With a luminosity of 1000 L Q (|Crimier et al.||2009[ hereafter 



Kempen et al. 2010| l. Sub-millimeter molecular line emission 



is a versatile probe of the physical and chemical conditions in 
these heated regions, the hot cores, revealing the initial con- 
ditions for forming stars and planetary systems. Hot cores are 
compact (sizes < O.lpc), warm (T > 100K), and show evidence 
of complex chemistry ( |Kurtz et al.|2000||Ceccarelli et al.|2007| l. 



* Herschel is an ESA space observatory with science instruments 
provided by European-led principal Investigator consortia and with im- 
portant participation from NASA 



Crim09) and a distance of only 440 pc ( |Hirota et al.||2007) > 
the protostar OMC2-FIR4 is an excellent laboratory to study 
hot core chemistry in the intermediate mass regime. A structure 
model, constrained by 7.5" to 14.8" resolution dust-continuum 
maps and the broadband spectral energy distribution, was made 
for OMC2-FIR4 by Crim09. The uncertainties in the proper- 
ties of the central component are large, but the model sug- 
gests a hot core radius of Ra-im = 440 AU and densities above 
~ 5 ■ 10 6 cm -3 . At ~1" resolution, FIR4 is seen to consist of 
several clumps, which may be forming a cluster of protostars 
( |Shimajiri et al.|2008 hereafter Shim08). Some of this activity 
may be triggered by an outflow from the nearby source OMC2- 
FIR3. 



SIMBAD entry: [MWZ90] OMC-2 FIR 4 
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Fig. 1. Averaged CH 3 OH line profiles in the E u < 200K (left) 
and E u > 200K (right) ranges, for bands lb (top), 2b (middle) 
and 3b (bottom). The gray area marks the biased weighted stan- 
dard deviation of the sample of normalized profiles. The two 
fitted Gaussians (blue, green) and their sum (dashed red line) are 
given. 



Methanol is a powerful diagnostic of the physical and chem- 
ical conditions in protost ellar sources (|van der Tak et aL] 2000 



|Leurini et alpOO?] [20071 |Maret et al.|2005||Wang et al.|2010) , 
and is used to that end here. A careful analysis of methanol is 
also important for recovering the other species which CH3OH 
lines often blend with. 



2. Observations 

The source OMC2-FIR4 was observed with the HIFI spec- 
trograph on the ESA Herschel Space Observatory (Pilbratt 
et al.|2010), u sing only the wide band spectrometer (WBS) ( de 
Graauw et alpO lO), covering bands lb (554.5 to 636.5 GHz), 
2b (714.1 to 801 GHz), and 3b (858.0 to 961.0 GHz). The WBS 
native resolution is 1.1MHz. The data were taken in dual beam 
switch (DBS) spectral scan mode with a redundancy of four. 
More observations, expected later in 2010, will cover most of 
HIFI's frequency range. Median baseline subtraction and side- 
band deconvolution were performed with the HIPE 2 software 
( |Ott|20 10). A frequency-dependent main beam efficiency within 
three percent of 0.70 and a calibration uncertainty of 15% were 
used in calculating line fluxes. For a detailed review of HIFI's 



in-orbit performance, see Roelfsema et al. (20101. 



The data quality in lb is superb, whereas 2b and 3b are sig- 
nificantly influenced by spurious features, "spurs", introduced 
by the local oscillator electronics, which act to increase noise 
in the single-sideband spectra and can damage spectral line pro- 
files. Combining the H and V polarizations, the RMS outside the 
bright lines is T* = 10 mK in lb, 21 mK in 2b, and 26 mK in 
3b. 



E„ range 


Fit 


FWHM [km/s] (bands lb, 2b, 3b) 


< 200 K 


1 

2 


(1.52, 1.84, 1.72) ± (0.01, 0.08, 0.08) 
(4.84, 5.16, 5.62) ± (0.03, 0.11, 0.07) 


> 200 K 


1 

2 


(2.22, 2.08, 2.06) ± (0.04, 0.11, 0.28) 
(5.68, 5.16, 4.44) ± (0.13, 0.13, 0.24) 




E„ range 


Fit 


Gaussian centre [km/s] 


< 200 K 


1 

2 


(0.02, -0.02, -0.06) ± (0.01, 0.03, 0.04) 
(-0.04, 0.03, -0.38) ± (0.01, 0.04, 0.04) 


> 200 K 


1 

2 


(0.05, 0.05, -0.22) ± (0.02, 0.04, 0.10) 
(0.39, 0.23, 0.04) ± (0.05, 0.05, 0.10) 



Table 1. Widths and central velocities of the two Gaussian fit 
components ("Fit") to each average line profile in Fig.[T] 



3. Results 

Preliminary data processing reveals 91 lines from 17 species in 
band lb, and similar numbers in 2b and 3b, establishing OMC2- 
FIR4 as a relatively line-rich protostar (see also the review by 
Ceccarelli et al. in this volume). Methanol lines are advanta- 
geous to study in OMC2-FIR4 because they are abundant, cover 
a large range of excitation conditions, and are straightforward to 
measure because the spectrum is less crowded than in higher lu- 
minosity hot sources, in which more levels of more species can 
be significantly excited. Out of the hundreds of CH3OH tran- 
sitions in each band, we detect 46 lines in band lb, 46 in 2b, 
and 40 in 3b, with a range of 40 K < E u < 540 K in excita- 
tion energy. Identification made use of the JPL databas^] with 
methanol data by |Xu et al. (2008). This paper discusses lines 
detected with > 4cr confidence and not blended or damaged by 
spurs. 



3. 1 . Line profiles 

Average line profiles were calculated in two energy ranges, 
E u < 200 K and E u > 200 K, for each observed band. These 
six average profiles are presented in Fig. [T] with the number of 
lines and relevant beam sizes given. Before averaging, each line 
profile was resampled to an 0.2 km/s resolution grid with origin 
at the rest frequency, fitted by a single Gaussian and normalized 
by the intensity at the fit center. Distinguishing between the core 
(|v| < 2 km/s) and wings ([v| > 2 km/s) of a profile, we note that 
the cores of the average line profiles seem to show a broaden- 
ing trend with decreasing beam size and increasing E u . Also, the 
wings are weaker in the higher energy regime than in the lower. 

Two-component blind Gaussian fits were made for the aver- 
age profiles in each energy range and band. The fits systemat- 
ically separate into a narrow and broad component, the widths 
and central locations of which are given in Table [T] The superi- 
ority of a two-Gaussian fit to that of a single Gaussian, as well as 
the trends mentioned in the previous paragraph, strongly suggest 
the presence of at least two physical emission components with 
different beam filling factors, as discussed in Sect. [4] 

3.2. Population diagram and LIE modelling 

The detection of 132 CH3OH lines covering 500K in upper level 
energy presents an excellent opportunity to perform a population 
diagram analysis (Goldsmith & Lange r|T999[ l to study the exci- 
tation conditions. In Fig. 2| we present the population diagram 



The catalogue is available on-line at http://spec.jpl.nasa.gov 
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Fig. 2. Observed and modelled (inset) CH3OH population dia- 
grams, plotting column density, divided by statistical weight, 
versus upper level energy for detected lines from all observed 
HIFI bands. Most transitions are given as crosses, but the K u = 
ladders for A and E states are indicated by filled squares and di- 
amonds, and K u = -3 by filled circles. The inset gives the syn- 
thetic population diagram from our LTE model for K u = and 
-3, with solid and open symbols representing the hot core and 
envelope contributions. The linear fit to the observed K u = A 
state lines is shown on both plots. 



for the observed methanol lines. Line fluxes were obtained from 
single-Gaussian fits. 

Part of the vertical scatter in Fig. [2] is due to optical depth, 
estimated from LTE modelling to contribute Aln(N u /g u ) < 
0.37. Another factor is beam dilution, which contributes up to 
Aln(N u /g u ) w 21n(37"/25") w 0.8 for a point source observed 
at the spectral edges of bands lb and 3b. However, most of 
the scatter is due to differences in thermalization conditions for 
CH3OH transitions with different upper level K quantum num- 
bers. We assume the K-ladders with the lowest n cl , e.g. K u = 0, 
approach LTE, as they thermalize below the densities expected 
for the hot core from the Crim09 model and are seen to lie close 
to each other, forming the upper envelope of the diagram. The 
K u = -3 ladder, which has a critical density of n cl as 10 9 cm -3 , 
forms the lower envelope, well below most detected ladders 
(-4 <K U < 4). 

Excitation temperatures and column densities from linear fit- 
ting are given in Table [2] for the ladder likely closest to LTE, 
K u - 0. The fit to the A state is shown in Fig. [2] and yields 
T rot = (145+12)Kand N A = (2.2 + 0.5) • 10 14 cnT 2 . Fitting the E 
state gives T rot = (120±8)Kand N A = (1.4 ± 0.3) • 10 I4 cnr 2 . A 
fit to the K u = -3 ladder gives T 10t = (66 ± 2)K. This low value 
is likely mostly due to subthermal excitation, but optical depth 
and beam dilution may also play a role, as explained earlier. 

A two-component LTE model was made with the CASSIS 
softwar^] treating methanol as a single species. The model en- 
velope and hot core have methanol abundances of 2 • 10~ 9 and 
10~ 6 , line widths of 3 km/s and apparent diameters of 20.0" and 
3.5", respectively. Source sizes, H2 column densities, and exci- 
tation temperatures were taken to be consistent with the Crim09 
source model, as well as with the observed T rot values in Table[2] 



Approach 


Component 


T TK1 (tVDe) 


N «,lcnr 2 l 

L MTieth L*- 111 J 


Population 


Kn = 0(A) 


145+12 (rot) 


(2.2 ± 0.5) • 10 14 


diagram 


K u = 0(E) 


120±8 (rot) 


(1.4 ±0.3)- 10 14 


LTE model 


Envelope 


40 (exc) 


2 • 10 14 


input 


Hot core 


120 (exc) 


6- 10 16 


LTE simu- 


K u = 0(A) 


164 (rot) 


1.4 • 10 14 


lated obs. 


K u = 0(E) 


168 (rot) 


1.7 ■ 10 14 



Table 2. From the top: the rotational temperatures and methanol 
column densities derived from the observed population diagram; 
the input excitation temperatures and column densities of the 
LTE model; the rotational temperatures and column densities de- 
rived from the modelled population diagram. 



and are discussed in Sect. [4] Table |2]gives the column densities 
and excitation temperatures we used as input to the model and 
the simulated observables derived from it. The simulated popu- 
lation diagram, including the Herschel-HIFI beam efficiency, is 
shown as the inset in Fig. [2] Optical depth in the modelled lines 
is typically 0.01 . . . 0.80. The LTE conditions align transitions 
from all /f-ladders, decreasing scatter in the simulated diagram. 

The model was adjusted to reproduce the K u = lines, which 
led to a solution where the hot core dominates the line fluxes, as 
seen in the inset of Fig. [2] As the T rot values in Table [2] show, 
a linear fit to the simulated observations yields results different 
from the model input. This is due to optical depth, blending, and 
beam dilution effects. The T rot derived from the simulation is 
sensitive to factor two changes in input line width at the 10% 
level due to optical depth effects. Owing to the high upper level 
energies of the detected K u = transitions, the envelope emis- 
sion only influences the derived T rot at the 1% level in the model. 
It is not clear that LTE is a valid assumption even for the K u = 
ladder, and thus a more comprehensive, non-LTE modelling ef- 
fort is underway. Preliminary results suggest the hot core density 
is higher than in the Crim09 best-fit model. 



4. Discussion 

The source OMC2-FIR4 is a line-rich protostar with a high 
degree of chemical and physical complexity. The hundreds of 
CH3OH transitions that are detectable with HIFI cover a wide 
range of excitation conditions and offer unprecedented spectral 
constraints on the structure of protostellar cores. 

Shim08, 



Previous knowledge (J0rgensen et al. 2006 



3 CASSI S has been developed by CESR-UPS/CNRS 
(http://cassis.cesr.frl. 



Crim09) suggests OMC2-FIR4 is dominated by three compo- 
nents: a large-scale, cool envelope ~ 10 4 AU across; a compact 
hot core ~ 10 2 AU across; and the outflow from FIR3, in par- 
ticular a suspected blue-shifted spot resulting from its interac- 
tion with the FIR4 envelope. The hot core gives rise to a domi- 
nant part of the emission in all lines in our LTE model, but the 
changes in line profiles through Fig.[T]imply reality is more com- 
plex. The envelope may contribute significantly to the lowest ex- 
citation lines, and shocks from regions such as where the FIR3 
outflow strikes FIR4 (Shim08) are likely to be important. 

The effects of the hot core radius and H2 column density on 
the line fluxes are degenerate if the lines are optically thin, as 
they are in our LTE model. Thus, an upper limit on the former 
and lower limit on the latter will result in an upper limit on the 
methanol abundance. Attributing the luminosity of FIR4, 10 3 L Q , 
to the central protostar and assuming grain mantle evaporation 
at 100 K, one finds an upper limit of R core ~ 760 AU for OA/im 
olivine grains. More realistic radiative transfer would decrease 
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this value. At 440 pc, 760 AU extends ~ 1.75" on the sky, which 
is adopted as the core radius in the model. 

To reproduce the observed K u = line fluxes with LTE 
at an assumed T exc = 120K, we need a column density of 
N(CH 3 OH) = 6 • 10 16 cirT 2 in the model hot core. Integrating 
the Crim09 H 2 density distribution from 100 to R Cr i m =440 AU, 
one obtains N(H2) ~ 6 ■ 10 22 cirT 2 . The source model was based 
on low-resolution maps, and the core is poorly constrained and 
excludes the inner 100 AU. We thus take the H2 column density 
as a lower limit and, adopting it in the LTE model, conclude that 
the obtained CH3OH abundance X core = 10~ 6 is a conservative 
upper li mit. The highest abund ance seen in low-mass hot cores is 
2 • 10~ 6 ( |j0rgensen et al.|2005 1. If the hot core does not dominate 
the line emission in reality, smaller hot core sizes and methanol 
column densities, as well as higher rotational temperatures, may 
be consistent with the data. Presently, the abundance of methanol 
in the envelope is very poorly constrained. We set it to 2 • 10" 9 in 
the LTE model, consistent with the factor 10 2 " 3 abundance jump 
seen in hot cores ( |van der Tak et al.|2000l|Maret et al.|2005) . 

The hot core density is broadly constrained to be between 
10 6 and 10 8 cirT 3 by this early analysis. Lower and higher val- 
ues would be difficult to reconcile with the large-scale density 
profile, furthermore the population diagram suggests the K lad- 
ders with the lowest n cr may approach LTE while those with the 
highest n cl clearly deviate from it. Typical kinetic temperatures 
well in the > 100 K hot core regime are implied by the obtained 
rotational temperatures of T rot w (120 . . . 145) K. Shock contri- 
butions from the FIR3-FIR4 interaction spot as well as the hot 
core itself need to be considered before drawing further conclu- 
sions about the temperature structure. 

An intriguing feature is the appearance of an asymmetry in 
the E u > 200 K lines, as seen in Fig. [T] Because these lines 
should be dominated by the hot core emission, the profiles ob- 
served could represent the blue asymmetric profile expected for 
infalling gas (e.g. Walker et al.|1994] |, implying we may be see- 
ing collapse in the central regions. To clarify whether infall, 
shocks, or other phenomena are responsible for the asymmetry, 
interferometry as well as further analysis of the Herschel-HIFI 
data, in particular the 13 CH30H lines, will be employed. 

We note the difference in column density of the K u — tran- 
sitions of the A and E states in Table|2] Whether it is truly signif- 
icant will be explored in a future paper. Also, early results from 
an analysis of formaldehyde in OMC2-FIR4 will be presented in 
a companion paper (Crimier et al., in prep.). 
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